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Excitation energies and transition probabilities of the first 2+ excitations in even tin and lead 
isotopes as well as the quadrupole moments of odd neighbors of these isotopes are calculated within 
the self-consistent Theory of Finite Fermi Systems based on the Energy Density Functional by 
Fayans et al. The effect of the density dependence of the effective pairing interaction is analyzed 
in detail by comparing results obtained with volume and surface pairing. The effect is found to be 
noticeable. For example, the 2"'"-energies are systematically higher at 200-300 keV for the volume 
paring as compared with the surface pairing case. But on the average both models reasonably agree 
with the data. Quadrupole moments of odd-neutron nuclei are very sensitive to the single-particle 
energy of the state A under consideration due to the Bogolyubov factor {u\ — v\). A reasonable 
agreement with experiment for the quadrupole moments has been obtained for the most part of odd 
nuclei considered. The method used gives a reliable possibility to predict quadrupole moments of 
unstable odd nuclei including very neutron rich ones. 

PACS numbers: 21.10.-k, 21.10.Jx, 21.10.Re, 21.60-n 



I. INTRODUCTION 

Presently there are two theoretical approaches which 
can quantitatively describe the bulk properties of nuclear 
isotope chains with a small number of effective coupling 
constants: selfconsistent mean field theories and density- 
functional theory. The successes and open problems of 
the mean field approaches are reviewed in Refs. [ll-Q- 
The Kohn-Shani density functional theory was originally 
proposed for chemistry and solids [J, Q • Important theo- 
retical developments have been made: an extension of the 
Hohenberg-Kohn theorem to pairing degrees of freedom 
by Oliveira, Gross, and Kohn allowed studies of super- 
fluids [li] and the generalization of functional theory to 
study excited states made it possible to investigate the 
electromagnetic response of correlated electron materials 
Q- In nuclear physics, a self-consistent Theory of Fi- 
nite Fermi Systems (TFFS) was derived by Khodel and 
Saperstein [8] on the basis of the TFFS by Migdal 9] sup- 
plemented with the many-body theory self-consistency 
relation llOll for the nucleon mass operator. As it was 
shown in [ll|, the self-consistent TFFS for nuclei with- 
out pairing can be reformulated as a particular version 
of the density functional method with a rather compli- 
cated p-dependence of the energy functional. It contains 
also r-dependent terms but with rather small strength 
resulting for the effective mass in a small difference of 
|to* „(r) — m\ ~ 0.05m. In a series of articles by Fayans 



et al. [il,[i3> (see also [Tj| and Refs. therein) the energy 
density functional (EDF) method was generalized for su- 
perfluid nuclei. Just as in the original Kohn-Sham ap- 
proach, the identity m* = m was imposed. A fractional 
form of the density dependence for the central part of 
the normal component of the EDF was introduced^ The 
coordinate dependence of it resembled that of 



111 but 



the functional form was much simpler making the self- 
consistent QRPA calculations easier. Note that a recent 
generalization of the Skyrme force in J15i | contains a new 
term with a density dependence resembling that in [ij]. 
In addition, the velocity dependent force in 15[ is rather 
weak leading to the effective mass m* ~ 0.9m. Thus, 
the selfconsistent mean field approaches may eventually 
converge with the density functional methods. 

The non-relativistic versions of the self-consistent 
mean field theories introduce three-nucleon forces which 
are often expressed as a density dependent two-body in- 
teraction. In general, one assumes a fractional power of 
the density dependence. Recent advances in Effective 
Field Theory open the possibility to connect the density 
functional with the effective two- and three- nucleon sys- 
tems which are determined from two-nucleon scattering 
and few-nucleon reactions. Reviews about the current 
status of such attempts are given in Refs.|16l.[l7t. 

The question arises whether the pairing interaction 
should have an analogous dependence on the normal nu- 
clear density. Several studies derived pairing interactions 



from free two-nucleon interactions. Baldo et al. solved 
the gap equation in semi- infinite nuclear matter [l8| , nu- 
clear slab [131 , and finite nuclei [20, [2l[ . The Paris and 
Argonne vig NN potentials were used, results being al- 
most identical. To make results more appropriate for 
practical nuclear self-consistent calculations dealing with 
pairing in a model space, the pairing problem was treated 
in a two-step way. The gap equation was solved in a 
model space with limiting energy Eq = 30 -^ 40 MeV 
with the use of the effective pairing interaction. The lat- 
ter is found in the subsidiary sub-space in terms of a free 
NN potential. For all systems under consideration and 
the two NN potentials the effective pairing interaction 
found is much stronger, up to ten times, at the surface 
than inside. The Milan group concentrated on the ^^^Sn 
nucleus, a traditional benchmark for the nuclear pairing 
problem, and solved the gap equation starting from the 
Argonne V14 potential ^^. In addition to the free NN 
interaction, they included corrections due to exchange 
with low- lying surface vibrations ("phonons") [2j] and 
high- lying excitations, mainly spin-dependent ones, [2J]. 
In the last article, a local 3-parameter density-dependent 
effective pairing interaction is constructed for the model 
space with Eq = 60 MeV which reproduce approximately 
exact gap values. Qualitatively, it is similar to that de- 
scribed above. Without all corrections, it consists of a 
strong surface attraction and very weak attraction inside. 
Taking into account of the phonon exchange makes the 
inner interaction repulsive. At last, inclusion of the spin- 
dependent excitations makes the inner repulsion rather 
strong. Thus, the ab initio calculations of the effective 
pairing interaction predict essential density dependence 
with strong surface attraction. 

As an alternative to consideration of the gap equation 
with complete realistic NN interactions, Bulgac and Yu 
used the fact that this equation depends mainly on the 
low-fc behavior of NN force which can be approximated 
with a rather simple analytical function. It helped to 
develop a renormalization scheme for the gap equation 
without any cutoff in terms of zero-range interactions 
with explicit coordinate dependence of the effective pair- 
ing interaction and to suggest an EDF for superfluid nu- 
clei [Uli^. 

The calculations by Fayans et al. employed both vol- 
ume pairing and surface pairing interactions. The bind- 
ing energies and the proton and neutron separation en- 
ergies were found to be insensitive to the type of pairing 
force used. But the odd-even staggering of charge radii 
can be quantitatively reproduced only if the strong den- 
sity dependence of the pairing force is introduced |14| . 

In this work, we investigate the excitation energies and 
transition rates of the low- lying 2"'"-states in spherical nu- 
clei with the aim to analyze the sensitivity of those ob- 
servables to the details of the pairing interaction. We will 
compare two opposite limits, the "volume pairing" with 
density independent effective pairing interaction J-^ and 
the case of the function JF^ with the surface dominance. 
The latter will be named for brevity the "surface pair- 



ing" . Several sets of calculations of these characteristics 
of the first 2"'"-excitations were carried out recently within 
the QRPA method with Skyrme force [23, ^^ and within 
the Generator Coordinate Method with the Gogny force 
[23 ] . No systematic analysis of the density dependence of 
the pairing force was performed in these studies. Deal- 
ing with low-laying quadrupole excitations, it is natural 
to include into analysis also quadrupole moments of odd 
nuclei which give test of static quadrupole pola rization. 
In this paper, we use the EDF method [lj| with the 
functional DF3-a [30|. In the latter the spin-orbit and 
effective tensor terms of the original functional DF3 [Ij, 
[ill were modified. All the QRPA-like TFFS equations 
are solved in the self-consistent basis {ex,(p\) obtained 
within the EDF method with the functional DF3-a. 



II. BRIEF OUTLINE OF THE FORMALISM 

For completeness, we describe shortly the EDF method 
of [lj| using mainly the notation of [3l[. In this method, 
the ground state energy of a nucleus is considered as a 
functional of normal and anomalous densities, 
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£[Pnir), pp{r),iyn{r),Vp{r)]d^r. 



(1) 



The normal part of the EDF fnorm contains the cen- 
tral, spin-orbit and effective tensor nuclear terms and 
Coulomb interaction term for protons. The main, cen- 
tral force, term of f norm is finite range with Yukawa-type 
coordinate dependence. It is convenient to extract the 
(5(r — r')-term from the Yukawa function separating the 
rest of 
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r — r 



-(5(r-r') (2) 



to generate the "surface" part f which vanishes in in- 
finite matter with p(r) = const. The Yukawa radius Vc 
is taken the same for the isoscalar and isovector chan- 
nels. The "volume" part of the EDF, £^{p), is taken 
in [13, HJ, im as a fractional function of densities p+ = 
Pn + Pp and p- = pn- Pp. 
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Here x — p^/{2po) is the dimensionless nuclear den- 
sity where po is the density of nucleons of one kind 
in equilibrium symmetric nuclear matter. The factor 
Co = {dn/deYy^ in Eq. ^ is the usual TFFS nor- 
malization factor, inverse density of states at the Fermi 
surface. 



To write down the surface term in a compact form 
similar to ^ , the "tilde" operator was introduced in [31[ 
denoting the following folding procedure: 
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D{r-r')(f>{r')dr'. 
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Then we obtain 



where 
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All the above parameters, aj.,a%, h\^ h\^ , h% , are dimen- 
sionless. 

In the momentum space, the operator ^ reads 



Dil) 



{qrc 



1 + {qr^y 



(8) 



In the small r^ limit it reduces to D{q) = — ((/t'c)^, and 
Eq. ([5]) could be simplified to a Skyrmc-likc form pro- 
portional to (Vp)^. 

The spin-orbit interaction reads 

J^sl = Corl{K+K'TlT2) [Vl(5(ri - r2) X (pi - P2)]-(cri+(T2), 

(9) 
where the factor rg is introduced to make the spin- 
orbit parameters k, k' dimensionless. It can be ex- 
pressed in terms of the above equilibrium density, Tq — 

In nuclei with partially occupied spin-orbit doublets, 
the so-called spin-orbit density exists. 



pliir) = Y,nl{cPl*ir){cTlWxir)), 



(10) 



where t = n,p — is the isotopic index and averaging 
over spin variables is carried out. As it is well known, see 
e.g. [8|, a new term appears in the spin-orbit mean field 
induced by the tensor forces and the first harmonic (fi of 
the spin Landau-Migdal (LM) amplitude. We combine 
those contributions into an effective tensor force or first 
spin harmonic, 

•^1 = Corligi + g[TiT2)S{ri - r2)((Ticr2)(piP2). (11) 

In Table 1, we present all parameters of the normal 
part of the EDF DF3-a we use. Note that the major 
part of these parameters is identical to the ones used in 
the DF3 functional [ij] . With one exception, all parame- 
ters for the central force part remained the same and only 
the spin-orbit and the first spin harmonic are changed ac- 
cording [201 • Application of the volume part dSj to equi- 
librium nuclear matter, with the equilibrium relation, i.e. 
vanishing pressure p{p) — p^d{E / p)/dp, permits to find 
the parameters a^ , h\^ and h\j^ in terms of the miclear 



matter density po, the chemical potential /io, and the 
compression modulus Kq = 9dp/dp. The asymmetry en- 
ergy parameter /3o yields a relation between the constants 
a^ , h^_ and /ij- • They are given in the upper half of Ta- 
ble 1. The radius vq introduced above is shown instead 
of Po- The value used in Ref. [iJ] was recalibrated in 
Ref. [30[ to obtain a more accurate description of nu- 
clear charge radii [32l | . One more remark should be made 
concerning the table. The "natural" TFFS normalizing 
factor Co = 2eoF/(3po) = 308.2 MeV fm'^ corresponding 
to parameters of nuclear matter in the third column of 
the tabic differs from the one, Cq — 300 MeV fm'^, recom- 
mended in the second edition of the Migdal's textbook on 
the TFFS [33] . To make a comparison with other articles 
within the TFFS, we recalculated all the strength param- 
eters to the latter. It explains a small difference of some 
values in the second column in the table from the original 
those in [iJ]. An essential difference between DF3 and 
DF3-a functional takes place for the "spin-dependent" 
sector in the bottom of the table. As we found in [30] , the 
second one describes the spin-orbit splitting of doublets 
better. 

The anomalous component of the EDF [14|reads 



£:a„(r) = ^-F«---(r;[p])l^-(r)l^ 
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where the effective pairing interaction reads: 

^« = Co/« = Co (.4 + h^^^^' + f^rliS/xf) . (13) 



The first two terms are similar to those in the TFFS 
^,^ or in the SHFB method jSQ]. The third in ^ is 
a new one introduced in [l3| . In this paper we use an ap- 
proximate version of (I13p with /y = 0. We will compare 
two models for nuclear pairing: the "volume" pairing 
{h^ = 0) and the "surface" pairing where both the pair- 
ing parameters /J^ and h^ are nonzero. One more remark 
should be made concerning the pairing problem. In the 
approach [l4| pairing was considered in the coordinate 
representation explicitly, solving the Gorkov equations 
with the method developed in Ref. [37|. However, it 
was found that the results are practically equivalent to 
those obtained within a more simple BCS-like scheme us- 
ing the gap A\\i = AxSw in a sufficiently large model 
space, e\ < -Emax- The effective pairing interaction (J13p 
for the BCS approximation is a little stronger than that 
in the coordinate representation (at ~ 5 -;- 10%, depend- 
ing on i?max)- For the systematic calculations in this 
article we use this simplified method of considering the 
pairing problem with -Bmax — 36 MeV. We do not apply 
this method for nuclei close to the dripline for which the 
diagonal approximation doesn't work [14]. 

Within the TFFS, the response of a nucleus to the ex- 
ternal quadrupole field Vq exp (iut) can be found in terms 
of the effective field. In systems with pairing correlations, 
equation for the effective field can be written in a com- 
pact form as 



Viiu)^VQ{iu)+fAioj)V{uj), 



(14) 



TABLE I: Parameters of the normal part of the EDF 



Parameter 


DF3 [14] 


DF3-a [30] 


MO, MeV 


-16.05 


-16.05 


ro, fm 


1.147 


1.145 


Ko, MeV 


200 


200 


/3, MeV 


28.7 


28.7 


a\ 


-6.598 


-6.575 


hl+ 


0.163 


0.163 


^2+ 


0.724 


0.725 


al 


5.565 


5.523 


hl- 








hl- 


3.0 


3.0 


a\ 


-11.4 


-11.1 


h% 


0.31 


0.31 


at 


-4.11 


-4.10 


ht 








Tc, fm 


0.35 


0.35 


K 


0.216 


0.190 


1 


0.077 


0.077 


91 








9. 


-0.123 


-0.308 



where all the terms are matrices. In the standard TFFS 
notation ,3, we have: 
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V^ = 







A(uS) = \ 0{uj) -Afiiuj) Af2{uj) 



(15) 



(16) 



(17) 



where /!, Mi, and so on stand for integrals over e of the 
products of different combinations of the Green function 
G(e) and two Gor'kov functios F^^) (e) and F^^) (e). They 
can be found in t9| and we write down here only the first 
of them which is of the main importance for us, 



C = 



de 
2Td 



G{e)G{e + u) - F^^^ (e)^^^) (e + u)] . (18) 



Isotopic indices in Eqs. (|15m7p are omitted for brevity. 
In Eq. (IT51) . F is the usual LM amplitude, 



F 



^ 

V 



(19) 



whereas the amplitudes F'^^ — F^" stand for the mixed 
second derivatives, 



J-"« 



SpSiy 



(20) 



In the case of volume pairing, we have F'^^ — 0. The 
explicit form of the above equations and (J18p is written 
down for the case of the electric (i-even) symmetry we 
deal with. A static moment of an odd nucleus can be 
found in terms of the diagonal matrix element (Ao|V(w = 
0)|Ao) of the effective field over the state Aq of the odd 
nucleon. 

The effective field operator ^(a;) has a pole in the ex- 
citation energy w^ of the state \s) under consideration. 



VQA{ujs)gQs] gos 
V{u;) = ^ ^ + Vr{oj). (21) 

U! — UJs 

The quantity g^s has the meaning of the corresponding 
excitation amplitude. It obeys the homogeneous coun- 
terpart of Eq. ([T4|) and is normalized as follows [3 , 



\+dA^ ' 



= -1, 



(22) 



with obvious notation. 

For excitation probabilities, it is more convenient to 
use the transition density operator which is conjugated 
to gos, 



Pos 



(23) 



The explicit definition of the normal and anomalous 
components of p}}'^ is as follows 



plf\r,r') = J^^6Gir,r';s,u,), 



tr(l,2)/ ,-, 
POs (r,!-) 



de 

2TTi 



,SF<^'''\r,r';e,u;s). 



(24) 



(25) 



The TFFS equation for transition densities for nuclei 
with pairing correlations. 



Pos = M^s)^Pas^ 



(26) 



is a complete analogue of the QRPA set of equations. 
Therefore we will often name it the QRPA equation. The 
transition density is normalized due to Eq. ([2^ . and the 
transition matrix element for the excitation of the state 
\s) with the external field Vq is given by 



Mos = J Vopl'Jr)dr. 



(27) 



III. 



CHARACTERISTICS OF THE 2+ 
EXCITATIONS 



The formalism described in the previous Section was 
used to describe 2f states in two isotopic chains of semi- 
magic nuclei, lead and tin. We investigate both a pure 
surface and a pure volume version of pairing. More cal- 
culational details can be found in Ref. J^]. We use 



the so-called developed pairing approximation. In par- 
ticular, we don't make any corrections to particle non- 
conservation effects induced with the Bogolyubov trans- 
formation. Therefore in the vicinity of double magic nu- 
clei, the results should be considered as very approxi- 
mate. As it was found in [ij], it is impossible to describe 
neutron and proton separation energies S'„ and Sp for all 
nuclei, from calcium up to lead, with sufficient accuracy 
using a fixed set of parameters in Eq. ((T3)) . the effective 
strength of the pairing interaction should be diminished 
with increasing nucleon number A. 
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FIG. 1: (Color online) Neutron separation energies Sn for 
lead isotopes. The volume pairing corresponds to (/|x ~ 
-0.31; h^ = 0), the surface one, to (/I, = -1.05; h^ = 0.94). 
The HFB theory predictions with the HFB-17 Skyrme func- 
tional are taken from ,38j . 



A exp. 
•- DF3-a, surf. 
o- DF3-a, vol. 
•- HFB-17 




100 105 110 115 120 125 130 135 
A 

FIG. 2: (Color online) Neutron separation energies Sn for 
tin isotopes. The volume pairing corresponds to (/^x ~ 
-0.33; ft* = 0), the surface one, to (/|„ = -1.05; h^ = 0.92). 
The HFB theory predictions with the HFB-17 Skyrme func- 
tional are taken from [3q ]. 

In this paper, we limit ourselves to two long isotopic 



chains, the lead and the tin chains. Therefore we deal 
with neutron pairing only. A short comment should be 
made on the procedure of solving the pairing problem. 
No particle number projection procedure is used in our 
calculations, i.e. particle number is conserved only on 
average, corresponding to the chosen chemical potential 
fi for the kind of nucleons under consideration. The accu- 
racy of this approximation is examined in a lot of papers. 
For the self-consistent SHE method with the SLy4 force, 
it was found in recent article [33| that the average differ- 
ence between exact and approximate gap values is 0.12 
MeV, the error being bigger in vicinity of magic nuclei. 

Eor finding the parameters of the pairing force (jI3p 
we use the strategy of "soft" variation of them to obtain 
better values of S'„ for both the chains under considera- 
tion. Values of Sn for both kinds of pairing are compared 
with the data in Eig. 1 and Eig. 2. Explicit values of 
the pairing parameters are given in the figure captions. 
Remind that we use the two-parameter version of (J13p . 
with /y = 0. For the volume pairing {h^ = 0), one pa- 
rameter remains which is smaller for lead than for tin 
approximately at 6%. Eor the surface pairing we deal 
with a two-parameter form of J-"^ . The "external" pair- 
ing parameter /l^ is taken ^-independent, in accordance 
with its physical meaning as the free NN T-matrix taken 
at negative energy £■ = 2/i flSJ. As to the second one, 
h^, it increases from the Sn chain to the Pb one at 2%, 
the resulting pairing attraction again becoming weaker, 
but only a little. Thus, the ^-independence of the pair- 
ing parameters in the case of surface pairing is weaker 
than for the volume one. This finding suggests to fa- 
vor surface pairing. As we see, the difference between 
the predictions for neutron separation energies is small 
for both versions and agreement with the experimental 
data is nearly perfect. Eor comparison, we display the 
pre dictions of the HEB-I7 version of the Skyrme force 
[36| which has a record accuracy in overall description 
of nuclear masses. We see that for these two chains our 
accuracy in description of neutron separation energies is 
even better. Of course, we achieved the agreement by a 
small variation of one of two paring parameters whereas 
calculations [3g are carried out with an universal set of 
parameters. However, the pairing part of the HEB-I7 
functional contains five parameters. 

Eig. 3 demonstrates that the normal neutron density 
/3„(r) and the anomalous one, t'n(r), both are practi- 
cally insensitive to the kind of paring used in the calcu- 
lation. On the contrary, the gap itself is very sensitive. 
For comparison, we took also a "medium" version, with 
(/Ix = —0.70; h^ = 0.50). It gives S'„ value approxi- 
mately with the same accuracy as the previous two. 

Let us now examine to what extent predictions for 
characteristics of 2f states are different for these two 
versions of pairing force which are equivalent in describ- 
ing the Sn values. A comment should be made be- 
fore presenting results of the QRPA calculations. Our 
QRPA code doesn't include the spin-orbit ^ and spin 
pri) terms of the effective interaction, therefore the self- 




r(fnn) 



FIG. 3; (Color online) Neutron density (pannel a), gap (b) 
and anomalous density (c) in ^""Pb nucleus. Solid and dot- 
ted lines correspond to the surface and volume versions cor- 
respondingly, the dashed one, to the medium version (/l^ = 
-0.70; /i« = 0.50). 



consistency is not complete and the excitation energy of 
the ghost l~-state does not automatically vanish. In the 
present investigation, we fine tuned the parameter a^ in 
([3]) in order to decouple the ghost state. This change is 
different for different nuclei but on average the value of 
|a^| increases at ~ 3% in comparison with that given in 
Table 1. 

Let us begin with the lead chain. Excitation ener- 
gies W2 are displayed in Fig. 4 and the probabilities 
B{E2,np), in Fig. 5. Experimental data for both quan- 
tities are taken from [40|. For comparison, results of 
the QRPA calculations of ^ with the SkM* and SLy4 
force are shown. Note that they were carried out with 
density independent pairing. We see that the difference 
Suj2 = w™' — w|"''^ is, on average, of the order of 0.3 
MeV, that is the effect under discussion is noticeable for 
this quantity. The results for volume pairing are sys- 
tematically higher, with the exception of the 2io,2i2p|-j 
isotopes for which the two versions practically coincide. 
Agreement with the data is, on average, quite reason- 
able for both the versions. Predictions of both the SkM* 
and SLy4 QRPA calculations for uj2 values have approx- 
imately the same accuracy as ours. 



For excitation probabilities the situation is more com- 
plex. For isotopes heavier than ^^*Pb our "surface" and 
"volume" curves are very close to each other. For lighter 
part of the chain the volume pairing generates larger 
probabilities than surface pairing does, producing differ- 
ences up to ~ 30%. Comparing with Fig. 4, we see 
that there is some unusual correlation between excita- 
tions energies and probabilities. Indeed, in magic nuclei 
where the pairing is absent for low-lying collective excita- 
tions there is a rule that a lower energy implies a larger 
probability. It can be qualitatively explained with the 
hydrodynamical Bohr-Mottelson (BM) model i4lj which 
gives a simple relation for the transition density of a L- 
vibration: 



tr,BM 

Pl 



dp 
ar 



(28) 



where ol = 1/^/2ujlBl, and Bl is the collective mass 
parameter of the BM model proportional to the nuclear 
mass. Then one obtains 



B{EL, np) 



2L + 1 
2ujlBl 



{Ml 



(29) 



where Mf^ 



{3Ze/4:TT)R'"^^, R being the nuclear ra- 
dius. Thus, in the BM model a lower value of the excita- 
tion energy wl inevitably leads to a higher value of the 
excitation probability. In our calculations, the situation 
is opposite. In principle, this is not strange. Indeed, even 
in magic nuclei the BM model works only qualitatively 
Q ■ If one solves equations of the self-consistent TFFS or 
any HF+RPA equations for nuclei without pairing, Eq. 
(|28p remains approximately true, but the mass parameter 
becomes w-dependent and deviates from the BM model 
prescription significantly '8i|. In nuclei with pairing, the 
situation becomes even more different from this simplest 
model as the normal component of the transition density 
(|24p depends now from the anomalous transition ampli- 
tudes cjQ^ ' (see Eq. (^51) ). They strongly depend on 
the kind of pairing. As a result, the correlation between 
the uiL and B{EL) values of the BM-type (|29l) can be 
destroyed. 

Experimental probabilities are known only for four 
even 204-2i0p^ isotopes. For all of them, the SkM* and 
SLy4 calculations are in perfect agreement with the data. 
Agreement of our calculations is poorer. It is especially 
true for the magic ^°^Pb nucleus where there is no any 
pairing. It should be noted that in this nucleus the collec- 
tivity of the 2^-state is not high: the B(E2) value is only 
about 8 single-particle units (spu). For a comparison, 
the B(E3) value for the 3]^-state exceeds 30 spu. But for 
excitations with low collectivity in nuclei without pairing 
the RPA solution depends strongly on the single-particle 
spectrum, and even a small inaccuracy in the positions 
of single-particle levels can change results significantly. 
In any case, some modification of the normal part of the 
functional DF3-a is necessary to obtain better agreement 
for the ^"®Pb nucleus. 



> 

CD 



LU 





' 


1 


1 — 


1 


' 


1 


1 1 1 1 1 1 1 


1 


1 1 1 


5 


- 














! 


- 




- 




Pb 






A exp. 


M 


- 




■ 












— •— surf. 


A 


■ 


4 


- 












— o— vol. 

--•--SkM* 

--♦--SLy4 


A 


- 


3 


~ 














/ 


: 


2 


- o^ 












i 


f 


/: 


1 


't~ 






=8;^ 


-^^ 


♦ ♦ t I— •— ♦"1 




fc-4: 




.▲ 


▲ 


▲ 


A 


▲ 


^*- 




▲ ▲ A. 




- 


1 


1 


1 




1 





1 


1 



184 



188 



192 



196 



200 



204 



208 212 



FIG. 4: (Color online) Excitation energies oj{2f) for lead isotopes. Predictions for mean field approach with the forces 
SkM*(dashed blue line) and SLy4(dashed green line) are taken from [23. The energy density functional results are given by 
the solid lines. 



In the tin chain, see Fig. 6, the situation with the ex- 
citation spectrum is partially similar to the one in lead. 
Again, 2^-levels are higher for volume pairing than in the 
surface case, and again the difference 6uj2 is — 300 keV, 
the surface predictions being closer to the experimental 
data. As to the SkM* spectrum, for isotopes heavier than 
^^^Sn it practically coincides with our "surface" one, both 
being higher than the experimental spectrum by approx- 
imately 200 -;- 300 keV. For lighter isotopes, it deviates 
from our surface spectrum significantly in an irregular 
way whereas the latter practically coincides with the ex- 
periment in this A region. As to the SLy4 spectrum, it 
also looks reasonable for the heavy part of the chain but 
for isotopes lighter of "'^^''^Sn it strongly oscillates around 
the experimental curve. In the dip minimum for ^^^Sn 
the UI2 value is less than the experimental one at approx- 
imately 1 MeV and it is close to an instability. 

The excitation probabilities are displayed in Fig. 7. 
Here the results show a very complex pattern. For the 
heavier part of the chain, beginning at the ^^"^Sn nucleus, 
our two theoretical curves and the SkM* practically co- 
incide, all being close to the experiment. The SkM* 
curve behaves in a non-regular way with strong devia- 
tions from the experimental data, up to ~ 50-^100%. The 
Sly4 interaction produces excitation probabilities which 
strongly decrease with the nucleon number A, implying 



drastic deviations from the data. The density functional 
approach is able to describe the A-dependence of the 
experimental B(E2,up) values rather well. For lighter 
tin isotopes, our two curves began to deviate from each 
other, the volume one being higher by ~ 25 -^ 30%, and a 
first glance may suggest that the volume pairing interac- 
tion performs much better. On the other hand, one has 
to notice the large error bars of the experimental data in 
the mass region below A=114. 

To investigate the role of pairing itself and of the type 
of its density dependence in detail, let us analyze dif- 
ferent components of the transition amplitude. Let us 
begin from the anomalous terms g^^'^-* (the index "Os" 
is for brevity omitted). They are displayed in Fig. 8 
for the ^™Pb nucleus. We see, first, that, for both the 
versions, the g^^' amplitude value is much bigger than 
Second, the coordinate dependence of the main 



(2) 



1^ . . 

g^^' amplitude is absolutely different for the two versions 
under comparison. In the surface pairing case, a strong 
surface maximum dominates whereas in the volume case 
g^^' is spread over the volume, with rather strong oscil- 
lations. In addition, it is seen that the integral effect of 
5surf should be noticeably bigger than that of g^J^. All 
this shows some asymmetry for Bogolyubov quasiparti- 
cles and quasiholes. Such a situation is typical for nuclei 
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FIG. 6: (Color online) Excitation energies aj(2^) for tin isotopes. Predictions for the SkM* and SLy4 force are taken from [27 
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which are close to the magic core. 




FIG. 8: (Color online) The neutron anomalous transition am- 
plitudes p'^'^' in ^""Pb nucleus. Solid lines correspond to sur- 
face, dotted to volume, and dashed, to the medium kind of 
pairing, see Fig. 3. 

The normal proton and neutron amplitudes g^^' for 
the same nucleus are displayed in Fig. 9. As we see, 
for this quantity the influence of the kind of pairing used 



is minimal. Thus, evidently, the rather big value of the 
difference 6102 — 300 keV for this nucleus is explained 
with different contributions of the anomalous amplitude 
g^^' which is much stronger in the case of surface pair- 
ing. For the transition densities, see Fig. 10, the effect is 
rather small but a little bigger than for the normal ampli- 
tudes g^^K This additional enhancement of the surface 
maximum of p*'"'-'^) (r) in the surface pairing case again 
originates from the term with g^^' in Eq. (I23p . In its 
turn, it explains the increase of the B{E2) value in this 
nucleus for the surface case. 

Let us go to the tin chain. Figs. 11-13 present for the 
^^®Sn nucleus the same quantities which were displayed 
in Figs. 8-10 for the ^oopj^ nucleus. This nucleus is in the 
middle of the chain, and all properties of the "developed" 
pairing, in particular, particle-hole symmetry should take 
place. Indeed, now (see Fig. 11) the amplitudes g'^-' 
and 5'^' possess a similar form and absolute value and, 
being of the opposite sign. In the result, we have \g^~' = 
g(i) _ 5(2)1 ^ y+) ^ g(i) + c,(2)| as it should be ^. 
Again, as in the 200p|-, case, the effect of the kind of 
pairing on the magnitude of g^^'^^ is drastic. As to that 
for the normal amplitudes g^^' and transition densities 



P 



fr{0) 



again it is rather moderate but of the another sign. 
Now in the volume case, the surface peaks in both these 
quantities are higher and, correspondingly, the B(E2) 
value is bigger. Evidently, in this case we deal with some 



10 








' ■ 












''/\ 


- 




^^^Sn 






■ 


- 




..g"" 


/ _, _ \ 


- 




"J^ 


- -_:. 




'^^^>~— tt™ 










.■■■■";^^^^^^, 






g"^' 




^ 








g'Vy 


■ 




■ 1 




1 i f i 1 i f 


i 1 i 



r{fm) 



r(fm) 



FIG. 9: (Color online) The proton and neutron normal transi- 
tion amplitudes g'"' in ^""Pb nucleus. Solid lines correspond 
to surface pairing, dotted ones, to volume pairing. 



FIG. 11: (Color online) The neutron anomalous transition 
amplitudes g^^'"^' in ^^*Sn nucleus. Solid lines correspond to 
surface pairing, dotted ones, to volume pairing. 
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FIG. 10: (Color online) The proton and neutron transition 
densities p*'"'''' in ^""Pb nucleus. Solid lines correspond to 
surface pairing, dotted ones, to volume pairing. 



FIG. 12: (Color online) The same as in Fig. 9 but for the 
"*Sn nucleus. 



destructive interference between normal and anomalous 
contributions to solutions of the equations of Section 2. 

To summarize, we see an effect of the type of pairing 
on the characteristics of the 2]^-states in spherical nuclei. 
The excitation energies l>ji are systematically lower in the 
surface case to ^ijJ2 — 300 keV, and the surface values are, 
as a rule, closer to the data. For B{E2) values, the effect 
is not so regular and here the volume version predictions 
on average look better. Thus, the present analysis is 
compatible with both volume and surface pairing. 

The issue could be naturally raised to what extent the 
small differences seen in the observables can be traced to 
the kind of pairing employed. In other words, is it possi- 
ble to fine tune the interaction parameters such that the 
both volume and surface pairing produce indistiguishable 
results, while keeping the mass differences and separation 



energies close to the experimental data? We carried out 
such an analysis for the tin chain. We consider the double 
mass differences 

D{N) - i (sn{N) - i (5„(iV - 1) + S,,{N + 1)) 

.(30) 
N even, which is very sensitive to the value of pairing 
gap. Note that the approximate relation D{N) ~ A 
takes place where A is an average gap value. We cal- 
culate the average difference between theoretical and ex- 
perimental values of this quantity, 



<^^) = jlJ E (^th W - i^cxp W)', (31) 

V N 



N even. We include into the analysis isotopes from '°Sn 
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FIG. 13: (Color online) The same as in Fig. 10 but for the 
"*Sn nucleus. 



FIG. 14; (Color online) The charge transition densities p*^^ ' 
in ^^*Sn nucleus. Solid lines correspond to surface pairing, 
dotted ones, to volume pairing. 



TABLE II: Dependence of the 2]''-state characteristics of the 
'^^^Sn nucleus on the strength of pairing force. 

version (SD) (MeV) L02 (MeV) B(£2, up)(e^b^) 



IV. QUADRUPOLE MOMENTS OF ODD 
NUCLEI 



surface 


0.10 


1.216 


0.172 


vol. f = -0.33 


0.11 


1.470 


0.206 


vol. f = -0.32 


0.16 


1.375 


0.193 


vol. f = -0.34 


0.11 


1.570 


0.216 



till "'^^^Sn for which the developed pairing approximation 
seems to be reasonable. Results are presented in Table 
II, for the surface pairing and for three versions of the 
volume pairing with different values of the strength pa- 
rameter /^. For all of them the characteristics of the 
2]^-state in the example ^-'^^Sn nucleus are given. It is 
seen that with increase of |/^| from the optimal value 
/^ = —0.33 deviations from the surface version predic- 
tions grow. With decrease of |/^| they become less, but 
this effect is much less than the initial difference even 
for the value /^ — —0.32 for which description of the 
mass differences is essentially worse than for the optimal 
value. An additional decrease of |/^| will absolutely de- 
stroy the nuclear mass description. In other isotopes of 
the tin chain, influence of variation of the /^ parameter 
to values of W2 and B{E2) is quite similar. Thus, the ef- 
fect under discussion originates mainly due to the surface 
nature of pairing versus the volume one. 

In conclusion of this Section we compare in Fig. 14 the 
charge transition density pli^{r) in the ^^^Sn nucleus with 
the experimental transition charge density found with a 
model independent analysis of the elastic electron scat- 
tering in [46|. The theoretical charge density is obtained 
from p*''(r) and p^ir) functions displayed in Fig. 13 with 
taking into account relativistic corrections [47[ . For both 
versions of pairing the agreement with the data is quite 
reasonable, and it is a little better in the surface case. 
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FIG. 15: (Color online) Static effective fields Vn Vp and d^ 
in ^"^Pb nucleus. Solid lines correspond to surface pairing, 
dotted ones, to volume pairing. 



Recently magnetic moments of odd spherical nuclei 
have been calculated [4^, |43] within the same self- 
consistent approach as the one used here. A reasonable 
description of the data for more than hundred of spher- 
ical nuclei was obtained. Especially high accuracy was 
reached for semi-magic nuclei considered in the "single- 
quasiparticle approximation" where one quasiparticle in 
the fixed state A = (n, l,j, m) with the energy e^ is added 
to the even-even core. According to the TFFS, a quasi- 
particle differs from a particle of the single-particle model 
in two respects. First, it possesses the local charge Cg (in 
our case, we have e^ = 1, e" = 0), and, second, the core 
is polarized due to the interaction between the particle 
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FIG. 16: (Color online) Static effective fields Vn Vp and d+ 
in ^^^Sn nucleus. Solid lines correspond to surface pairing, 
dotted ones, to volume pairing. 



and the core nucleons via the LM amplitude. In other 
words, the quasiparticle possesses the effective charge Ceff 
caused by the polarizability of the core which is found by 
solving the above TFFS equations. In the many-particle 
Shell Model [50|, a similar quantity is introduced as a 
phenomenological parameter which describes polarizabil- 
ity of the core consisting of outside nucleons. 

In non-magic nuclei, the term quasiparticle takes a 
double meaning. In addition to the initial LM concept 
we consider the Bogolyubov quasiparticles with occu- 
pation numbers nf={Ex — e\)/2E\ and energies E\ — 
\/{e\ — y^Y -\- A| and solve the set of the QRPA equa- 
tions (14) instead of one RPA equation. 

The success of the single-quasiparticle approximation 
in describing the magnetic moments of semi-magic nuclei 
makes it of interest to try to use the same approach for 
quadrupole moments. In this article, we do such anal- 
ysis limiting ourselves with odd neighbors of the even 
tin and lead isotopes considered in the previous Section. 
To our knowledge, there is no systematic calculations of 
quadrupole moments of these nuclei. 

The static quadrupole moment of an odd nucleus in 
the single particle state A can be found in terms of the 
effective field ([T4)) with the static external field Vb = 
^167r/5r2y2o(6') as follows ^^■. 



where ux, vx are the Bogolyubov coefficients and 



Vx 



2j 



Vir)Rii^{rydr. 



(32) 



(33) 



The j-dependent factor in ([55]) appears due to the angu- 
lar integral 52|. For j > 1/2 it is always negative. For 
odd neighbors of a magic nucleus the "Bogolyubov" fac- 
tor in ([5^ reduces to 1 for a particle state and to — 1 for 
a hole one. 



is V"'P{r) and d+(r)=d„ (r) -I- dn {r), are displayed in 
Figs. 15, 16 for ^°*Pb and ^^^Sn nuclei, correspondingly. 
Note that the identity d~ {u}=0)=0 takes place |9|. One 
can see large surface maxima of the quantities V"'P(r) 
similar to those in Figs. 9, 12 for the BM-like transition 
amplitudes gn/p{r). In-volume ("quantum") corrections 
are relatively small, therefore the integral in Eq. ([55)) 
is always positive. For protons, it is noticeably larger 
than the similar integral with the bare field V'^, see the 
discussion on the effective charges below. 

Diagonal matrix elements l\'6'3\i of the proton effective 
field are displayed in Fig. 17 for the tin isotopes and in 
Fig. 18, for the lead ones. As it is seen, for a major part of 
the tin isotopes, the difference between values of proton 
matrix elements V^ surface and volume pairing is quite 
small. Only for ^^^^^^^Sn nuclei it reaches 10%. In the 
lead region, the difference is more pronounced reaching 
~ 30 -=- 40% for 9/2" and 11/2" states. 

Corresponding quadrupole moments for nuclei with 
odd proton number Z = 50 ± 1 and Z = 82 ± 1 are 
presented in Table III. As it was noted above, in this 
case the Bogolyubov factor in ([5^ is trivial, equal to ±1. 
In order to check our approach, we selected only nuclei 
where there are experimental data and those which sat- 
isfy presumably the single-quasiparticle approximation. 
In particular, we excluded several light Tl isotopes with 
known quadrupole moments of low-lying excited 9/2~ 
states. If to suppose that they are single-quasiparticle 
l/ig/2 states, they should have essentially higher excita- 
tion energies than it takes place. 

Experimental data are taken from the compilation [53] . 
From several cases of proton excited isomeric states we 
limit ourselves with only two, the Ig^/j state in the "'^^^Sb 
and 2dJ^,2 state in ^os^p^ nuclei, for which the hypothe- 
sis on the single-quasiparticle structure seems to us more 
or less safe. Again, we presented results for both the 
kinds of nuclear pairing (the quantities Q^JJ'^ and Q^^' 
for surface and volume pairing, correspondingly). In the 
6-th column of the table, the single-particle quadrupole 
moment is presented which is found from Eqs. (|5^ . ([55)) 
with substitution V — >■ Vq. As it follows from Fig. 17, for 
odd- proton neighbors of the tin isotopes, difference be- 
tween values of quadrupole moments for surface and vol- 
ume pairing is quite small, in limits of 10%. In the lead 
region, see Fig. 18, the difference is more pronounced, 
but here the number of the data is very small, only 4. 
In addition, only in the 203,205gj g^j^^ 205^^^ ^g^gg ng^tron 

pairing exists. For these nuclei, the effect under discus- 
sion reaches ~ 30 -^- 40%. 

For the long chain of twelve In isotopes agreement with 
the data is quite reasonable. For five Sb isotopes (six val- 
ues of the quadrupole moment) agreement is rather poor, 
disagreement reaching ~ 50 -;- 100%. A similar situation 
takes place for two lighter Bi isotopes. For the ^"^Bi iso- 
tope where pairing is absent experimental data are con- 
tradictory. We think that the main reason of existing 
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TABLE III: 
clei. 



Quadrupole moments Q (e b) of odd-proton nu- 



nucl 



Ao 



?surt 
th 



^vol 



surt vol 
ficff fioff 



107 
109 
111 
113 
115 
117 
119 
121 
123 
125 

i"Sb 
"^Sb 

121 



139/2 
I59/2 
I59/2 
lff9/2 
lS9/2 
lff9/2 
I99/2 
lS9/2 
lff9/2 
lS9/2 
I99/2 
I59/2 



2d 

2d; 

2d 



Sb 2d 



■5/2 

/2 
■5/2 
'5/2 



'Sb 



l57*/2 
lS7/2 

3d. 



'Bi l/i. 



'Bi 
'Bi 



'3/2 
9/2 



l/i9/2 
I/I9/2 



-H0.83(5) 

-f0.81(5) 

-1-0.84(3) 

-f0.80(2) 

-1-0.80(4) 

+0.81(5), 0.58(9) 

+0.829(10) 

+0.854(7) 

+0.814(11) 

+0.757(9) 

+0.71(4) 

+0.59(3) 

-0.36(6) 

-0(2) 

-0.37(6) 

-0.36(4), -0.45(3) 

-0.48(5) 

-0.49(5) 

0.74(15) 

-0.68(6) 

-0.59(4) 

-0.37(3), -0.55(1) 

-0.77(1), -0.40(5) 



+0.83 
+0.98 
+1.11 
+1.16 
+1.12 
+1.03 
+0.96 
+0.91 
+0.83 
+0.74 
+0.66 
+0.55 
-0.88 
-0.82 
-0.77 
-0.72 
-0.81 
-0.74 
+0.23 
-1.32 
-0.94 
-0.34 



+0.90 
+1.07 
+1.14 
+1.10 
+1.02 
+0.97 
+0.95 
+0.92 
+0.84 
+0.74 
+0.74 
+0.49 
-0.81 
-0.77 
-0.76 
-0.73 
-0.81 
-0.74 
+0.23 
-0.91 
-0.72 
-0.34 



+0.18 
+0.18 
+0.18 
+0.19 
+0.19 
+0.19 
+0.19 
+0.19 
+0.19 
+0.19 
+0.19 
+0.19 
-0.14 
-0.14 
-0.14 
-0.14 
-0.17 
-0.17 
+0.12 
-0.25 
-0.25 
-0.25 



4.6 
5.4 
6.2 
6.1 
5.9 
5.4 
5.1 
4.8 
4.4 
3.9 
3.8 
2.9 
6.3 
5.9 
5.5 
5.1 
4.8 
4.4 
1.9 
5.3 
3.8 
1.4 



5.0 
5.9 
6.3 
5.8 
5.4 
5.1 
5.0 
4.8 
4.4 
3.9 
3.9 
2.6 
5.8 
5.5 
5.4 
5.2 
4.8 
4.4 
1.9 
3.6 
2.9 
1.4 



disagreements is neglecting the phonon coupling effects. 

Let us go to odd-neutron nuclei, the odd tin and lead 
isotopes. The results are presented in Table IV and Figs. 
19,20. In selecting nuclei for the table, we used the same 
concept as for protons. In this case, we included into the 
analysis twelve excited states, in addition to the ground 
ones. With the only exception of the ^"^Pb nucleus, all 
the nuclei under consideration exhibit pairing effects and 
the factor {u\ — v\) in Eq. ([5^ becomes non-trivial. It 
changes permanently depending on the state A and the 
nucleus under consideration. Note that in the case of 
magnetic moments the factor of (u? + '^^a) = ^ appears 



in the relation analogous to (|32p |51| . In our case, this 
factor determines the sign of the quadrupole moment. 
In all cases when the sign of the experimental moment is 
known the theoretical sign is correct. This permits to use 
our predictions to determine the sign when it is unknown. 
The factor under discussion depends essentially on values 
of the single-particle basis energies £\ reckoned from the 
chemical potential /z as we have {u\ — v\) = {e\ — 11) /E\. 
Keeping in mind such sensitivity, we found this quantity 
for a given odd nucleus {Z,N + 1), A^ even, with taking 
into account the blocking effect in the pairing problem 
51| putting the odd neutron to the state A under con- 
sideration. For the Vx value in Eq. (15^ we used the 
half-sum of these values in two neighboring even nuclei. 
We consider agreement with the data reasonable if we 
have IQth — Qexpl < 0.1 + 0.2 e b. If to use such a 
criterion, there are 7 "bad" cases in Table IV, and 16 
"good" . Several rather strong disagreements with the 



experimental data in Table IV for high-j levels l/iii/2 in 
Sn isotopes and IJ13/2 in Pb isotopes originate just from 
their too distant positions from the Fermi level. Thus, 
the Q values depend strongly on the single-particle level 
structure. Again, as for protons, the difference between 
predictions of the two models under consideration is, as 
a rule, rather small, and only for li]^3 /2-states in the lead 
chain it reaches ~ 20 + 30%. 

In the last two columns of Tables III and IV, the ef- 
fective charges are presented which are defined as e^^ = 
V^'^/{Vq)\. It is a direct characteristic of the core polar- 
izability by the quadrupole external field. In these tables, 
there are only two nuclei, ^"^Bi and ^°^Pb, with a double- 
magic core, and in this case the polarizability is rather 
moderate, e^^j = 1.4, e"fj = 0.9. In nuclei with an unfilled 
neutron shell it becomes much stronger, eoff ~ 3 + 6. 
The reason is rather obvious. Indeed, for the case of a 
positive parity field Vq, virtual transitions inside the un- 
filled shell begin to contribute and small energy denomi- 
nators appear in the propagator £" p8|) playing the main 
role in the problem under consideration. It enhances the 
neutron response to the field Vq and, via the strong LM 
neutron-proton interaction amplitude J^"p, the proton re- 
sponse as well. Results for the chain 203,205,209gj show 
how the polarizability grows with increase of the number 
of neutron holes. Keeping in mind this physics, one can 



represent the effective charges as e^j 



eff 



l+Cpoi' e. 



eff 



--pol 



where Cp^j is the pure polarizability charge. To separate 
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FIG. 17: (Color online) Diagonal matrix elements V^ of the effective proton quadrupole field in the tin isotopes. Solid lines 
correspond to surface pairing, dotted ones, to volume pairing. 
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FIG. 18: (Color online) The same as in Fig. 17, but for the lead isotopes. 
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FIG. 19: (Color online) Quadrupole moments of odd tin isotopes. Solid lines correspond to surface pairing, dotted ones, to 
volume pairing. Experimental data are shown with A for 3/2+, T for ll/2~, A for 5/2"'', and V for 7/2"*" states. 
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FIG. 20: (Color online) Quadrupole moments of odd lead isotopes. Solid lines correspond to surface pairing, dotted ones, to 
volume pairing. Experimental data are shown with A for 13/2+ , ▼ for 3/2", and A for 5/2~ states. 
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contributions of the unfilled shells and core nucleons ex- 
plicitly, one can divide the Hilbert space of the QRPA 
equations (fT4| to the "valent" and subsidiary ones and 
carry out the corresponding renormalization procedure 



V. DISCUSSION AND CONCLUSIONS 

The effect of the density dependence of the pairing 
interaction to low-lying quadrupole excitations in spher- 
ical nuclei is analyzed for two isotopic chains of semi- 
magic nuclei. Static quadrupole moments of neighbor- 
ing odd nuclei are also examined. The complete set of 
the QRPA-like TFFS equations for response functions is 
solved in a self-consistent way within the EDF approach 
to superfluid nuclei with previously fixed parameters of 
the functional. The DF3-a functional [301 i^ used which 
is a small modification of the functional DF3 [1^ [l4| . 
Specifically, spin-orbit and effective tensor terms of the 
initial EDF DF3 were changed. Two models for effective 
pairing force are considered, the surface and the volume 
ones, which give rise to approximately the same accu- 
racy in reproducing mass differences. A noticeable effect 
in excitation energies UJ2 is found: predictions for the 
volume model are systematically higher than the surface 
ones by 5u!2 — 200 -;- 300 keV. As to the excitation prob- 
abilities B{E2,-wp), the effect is not so regular, however, 
as a rule, the volume values are also higher. Thus, the 
correlation between these two quantities typical for the 
BM model, where a higher frequency always results in 
a lower probability, is destroyed. On the average, both 
models reasonably agree with the data. In addition, they 
both reproduce rather well the model-independent charge 
density /3'fj(2^) for the ^-^^Sn nucleus. 

Comparison with recent QRPA calculations [23] with 
the Skyrme force SkM* and SLy4 shows that for the lead 
chain they agree with the data a little better than our 
results but for the tin chain the situation is opposite and 
our predictions occur to be essentially better. The sur- 
face model is systematically better in describing the en- 
ergies W2 whereas the excitation probabilities are, as a 
rule, reproduced better with the volume model. 

Whereas the c harg e radii study [ij] and ab initio the- 
ory of paring [l3,[2Jl favor the surface pairing, the w(2^) 
and B{E2, up) data do not allow to prefer any of the two 
kinds of pairing. 

A reasonable agreement with experiment for the 
quadrupole moments of odd neighbors of the even tin 
and lead isotopes has been obtained for the most part 
of nuclei considered. For odd-proton this confirms that 
the single-quasiparticle approximation works sufficiently 
well. For odd-neutron isotopes under consideration, va- 
lidity of this approach was checked previously with the 
analysis of magnetic moments |49| . In the case we con- 
sider, the problem is more complicated than for odd 
proton isotopes as the Bogolyubov factor (u| — V\) = 
{e\ — ^)/E\ comes to the quadrupole moment value, in 



addition to the matrix element of the effective field V\. 
This factor makes the quadrupole moment value very sen- 
sitive to accuracy of calculating the single-particle energy 
e\ of the state under consideration, especially near the 
Fermi surface as the quantity Q\ vanishes at e\ = fi. For 
such a situation, the influence of the coupling of single- 
particle degrees of freedom with phonons, see [i^, [55j . 
should be especially important. This rather complicated 
problem will be considered separately. 

As to the effect of the density dependence of pairing, 
for quadrupole moments it is, on the average, less than 
for quadrupole transitions. It depends on a nucleus ex- 
amined and on the odd-nucleon state as well. In the tin 
region, it is, as a rule, of the order of ~ 10%. However, 
in the lead region it is higher and reaches ~ 30 -^ 50% for 
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